/xSR studies of the frustrated quasi-2d square-lattice spin system 
Cu(Cl,Br)La(Nb,Ta) 2 07: evolution from spin-gap to antiferromagnetic state 
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We report muon spin relaxation (/iSR) and magnetic susceptibility measurements on 
Cu(Cl,Br)La(Nb,Ta)207, which demonstrate: (a) the absence of static magnetism in 
(CuCl)LaNb207 down to 15 mK confirming a spin-gapped ground state; (b) phase separation be- 
tween partial volumes with a spin-gap and static magnetism in (CuCl)La(Nb,Ta)20r; (c) history- 
dependent magnetization in the (Nb,Ta) and (CI, Br) substitution systems; (d) a uniform long-range 
collinear antiferromagnetic state in (CuBr)LaNb20?; and (e) a decrease of Neel temperature with 
decreasing Br concentration x in Cu(Cli_ ;E Br ;E )LaNb207 with no change in the ordered Cu moment 
size for 0.33 < x < 1. Together with several other /xSR studies of quantum phase transitions in 
geometrically-frustrated spin systems, the present results reveal that the evolution from a spin-gap 
to a magnetically ordered state is often associated with phase separation and/or a first order phase 
transition. 

PACS numbers: 75.35.Kz 73.43.Nq 76.75.+i 



I. INTRODUCTION 

Modern studies of quantum phase transitions (QPTs) 
seek novel features of ground states near phase bound- 
aries. In a narrow range of J2/J1 ratios of square lat- 
tice spin systems, geometrical frustration of the nearest 
neighbour (Ji) and next nearest neighbour (J2) exchange 
interactions is expected to yield a spin-gap state. Recent 
synthesis of a new square lattice system (CuCl)LaNb2 07 
brought the first long-awaited system which might help 
elucidating this hypothesis and the associated QPT. In 
this paper, we report muon spin relaxation (/iSR) and 
low-field susceptibility studies of this system and rele- 
vant compounds obtained by (CI, Br) and (Nb,Ta) sub- 
stitutions. 

Unlike thermal phase transitions which represent 
changes of systems as a function of temperature, stud- 
ies of QPTs follow the evolution of ground states across 
phase boundaries at T — > by varying, for example, 
chemical composition or pressure as a tuning param- 
eter. Recent experimental studies revealed novel and 
sometimes unexpected features of QPTs, such as first- 
order- like behavior, phase separation and slow spin fluc- 
tuations at boundaries of magnetically-ordered and para- 
magnetic states in itinerant electron magnets MnSA^ 
and (Sr,Ca)Ru03 3 , which are metallic systems without 
spin frustration. Phase separation has also been found 
between the stripe spin-charge ordered state and the su- 
perconducting states without static magnetism in high- 
T c cuprate systems^*^ involving holes doped in anti- 
ferromagnetic Cu02 planes. It is interesting to com- 
pare these results to QPTs of insulating and geomet- 



rically frustrated spin systems (GFSS) 7,8 at boundaries 
between magnetically ordered states and spin-gap/spin- 
liquid states realized without static magnetism. 

GFSS on triangular or Kagome lattices exclusively in- 
volve antiferromagnetic nearest-neighbour interactions in 
hexagonal geometry. In contrast, the two-dimensional (2- 
d) square-lattice J1-J2 system is unique in its underlying 
square lattice geometry as well as the possible involve- 
ment of ferromagnetic interactions. Several compounds 
of vanadium oxide&i , synthesized in the early search 
for spin-gapped J1-J2 materials, unfortunately showed 
magnetic order at low temperatures (2.1-3.5 K). As il- 
lustrated in Fig. lc, the J2/J1 ratios of these systems, 
estimated from magnetic susceptibility, narrowly missed 
the region predicted for formation of a singlet ground 
state with a spin-ga p — 

In 2005, Kageyama and co-workers 11 succeeded 
in synthesizing a spin-gap candidate J1-J2 system 
(CuCl)LaNb207, which has a quasi 2-d crystal struc- 
ture with S=l/2 Cu moments as shown in Figs, la 
and b. The magnetic susceptibility \ of this system, 
shown in Fig. 2, exhibits typical spin-gap behavior with 
an estimated gap A/ks = 27 K. The spin gap was 
directly observed by inelastic neutron scattering 11 and 
high-field magnet iz at ioni 2 -. The very small Weiss tem- 
perature 6 = J\ + J2 < 5 K in the 1/x plot (Fig. 2) 
indicates that J\ and J2 have nearly equal magnitudes 
but opposite signs. Based on these observations, we as- 
sign this compound to fall in the spin-gap region in Fig. 
lc. 

A sister compound (CuBr)LaNb207 orders into a 
collinear antiferromagnet (CLAF) below T/v = 32 K, as 
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FIG. 1: (color) (a) Crystal structure of (CuCl)LaNb 2 07-^. 
(b) Exchange interactions J\ and J 2 on the 2-d square lat- 
tice, (c) Conceptual phase diagram of the spin 1/2 square 
lattice J1-J2 model 14,15 as a function of J\ and J 2 with re- 
gions of collinear antiferromagnetic (CLAF), ferromagnetic 
(FM) and Neel antiferromagnetic (NAF) order. A spin-gap 
or spin-liquid state is expected in the shaded region. The 
present work elucidates the evolution illustrated by the pur- 
ple arrow. 
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FIG. 2: (color) Magnetic susceptibility \ of (CuCl)LaNb 2 7 , 
(CuCl)LaTa2 07 and % and 1/% of (CuBr)LaNb207 in an ex- 
ternal field of 1 kG. 



shown by neutron scattering 13 and susceptibility mea- 
surements (Fig. 2). This indicates a slight change of the 
J2 1 J\ ratio of the Br compound from that of the CI com- 
pound (Fig. lc), together with a relatively large magni- 
tude of the antiferromagnetic J2 which is at least enough 
to support this Neel temperature. In order to study 
the evolution from the CLAF to the spin-gap states, il- 
lustrated by a purple arrow in Fig. lc, we investigate 
solid solution systems C^Cli-^Br^LaM^Oy as well as 
(CuCl)La(Nbi_ 2y Ta 2/ ) 2 ? 1 ^. The latter substitution is of 
great interest since it tends to suppress the spin gap, as 
seen in x for ?/ = 1 in Fig. 2, without perturbing the 
direct exchange path on the CuCl square lattice plane. 
This may allow finer tuning of the J2/J1 ratios with 
smaller effect of randomness. 

We started a project to study these systems by /iSR 
in 2005. Since then a few parallel efforts have been 



made in neutron scattering, high-field susceptibility / 
magnetization, NMR, and Raman scattering measure- 
ments. While NMR and Raman studies have been / 
will be reported in separate independent papersiii 8 /iSR 
and the remaining studies and characterization of the 
specimens will be reported in three consecutive papers 
(I, II 20 and III 19 ). In this paper (I), we focus on /iSR 
and low-field magnetic susceptibility measurements on 
Cu(Cl,Br)La(Nb,Ta)207, which demonstrate: (a) ab- 
sence of static magnetism in (CuCl)LaNb2 07 down to 15 
mK confirming the spin-gapped ground state; (b) phase 
separation between partial volumes with a spin-gap and 
static magnetism in (CuCl)La(Nb,Ta)2 0r; (c) history- 
dependent magnetization in the (Nb,Ta) and (CI, Br) sub- 
stitution systems; and (d) a uniform long-range collinear 
antiferromagnetic state in (CuBr)LaNb207. Details of 
preparation of specimens, elastic neutron scattering and 
high- field magnetization/susceptibility studies will be re- 
ported in paper II for (Nb,Ta) substitutions and III for 
(CI, Br) substitutions. The companion paper II is sub- 
mitted simultaneously with the present paper so that it 
may be published back- to-back with the present paper. 
Paper III will be published separately. 

The development of theories for J1-J2 systems are in 
progress. Following initial conjectur e - 14 ! 15 of existence of 
a spin gap in the border region of CLAF state and fer- 
romagnetic (FM) state, a more recent theory 21 predicts 
nematic spin arrangement, instead of the spin gap, at 
this border, when first and second neighbour interactions 
are considered in a Heisenberg model of a 2-dimensional 
square lattice. However, a spin-gap state could well be 
expected for situations involving higher order interac- 
tions, such as, 3rd and 4th nearest neighbour exchange 
interactions, anisotropy or three-dimensional interlayer 
couplings. In a recent experimental effort to develop com- 
pounds relevant to the present Cu(Cl,Br)La(Nb,Ta)207 
system, Tsujimoto et al^ synthesized (CuBr)A 2 B 3 Oio 
(A = Ca, Sr, Ba, Pb; B = Nb, Ta) which has three per- 
ovskite layers between the adjacent CuBr planes. One 
of these systems, (CuBr)Sr2Nb30io, exhibits a positive 
Curie temperature when the high-temperature suscepti- 
bility x is extrapolated to low temperatures in a plot of 
1/x versus T. This suggests that a fine tuning of parame- 
ters indeed brings the square lattice Cu(Cl,Br) plane into 
the side of FM correlations across the spin gap region, as 
illustrated in Fig. lc. 

In an NMR study of CuClLaNb 2 7 Yoshida et al£L 
found a signature of dimerization, although the corre- 
sponding superlattice satellite peaks of X-ray scattering 
are rather weak in intensity, and the assigned symme- 
try is not fully consistent with the results of recent Ra- 
man measurements 18 . In view of these developments, 
here we adopt the J1-J2 model as an appropriate start- 
ing point for discussions of Cu(Cl,Br)La(Nb,Ta)207, al- 
though there may be some influence of higher order in- 
teractions, dimerization, and other effects existing in real 
materials. Even when a small dimerization is essential 
for the formation of the spin gap, all the experimental 
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results described in this paper can still be regarded as 
elucidating an interesting boundary between a spin-gap 
and CLAF state. 

During the course of this study, we found a his- 
tory dependence of the magnetic susceptibility of 
Cu(Cl,Br)La(Nb,Ta) 2 7 which sets in at T - 7 K in 
a low applied field of ~ 100 G. This feature will be re- 
ported in section IV, following the /iSR results in section 
III. A renewed interest in J\- J 2 systems was generated by 
the recent discovery of La(F,0)FeAs superconductors 23 
which have Fe moments in this geometr y 24 ! 25 . In sec- 
tion V, we will compare the present results with fiSR 
studies of other GFSS, including Kagome lattice systems, 
cuprate and other systems on body-centered tetragonal 
lattices, and FeAs superconductors. 

II. EXPERIMENTAL METHODS 

Polycrystal specimens of solid solution systems 
Cu(Cli_ x Br x )LaNb 2 7 and (CuCl)La(Nbi_ 2y Ta 2y ) 2 7 i^ 
were synthesized at Kyoto University using ion-exchange 
reactions at low temperatures, as described in refs j 19 i 2Q . 
X-ray diffraction results show no trace of impurity phases 
for the entire concentration regions < x < 1 and 
< y < 1 within the experimental detection limit. In 
the case of (CI, Br) substitutions, the homogeneous and 
random distribution of substituted atoms has been ver- 
ified by the variation of the a- and c-axis lattice con- 
stants against composition x 19 following Vegard's law. 
For the (Nb,Ta) substitutions, similar information could 
not be obtained from X-ray diffraction due to nearly 
equal atomic radius of Nb and Ta. Chemical homogene- 
ity of both the (Br, CI) and (Nb,Ta) samples have been 
checked using energy dispersive spectroscopy (EDS) of 
transmission electron microscope (TEM) measurements, 
which confirmed uniform solutions with the spatial res- 
olution of 10 n m 19 ' 20 , The specimens were pressed into 
disc-shaped pellets having typical dimensions of 10 mm 
in diameter and 1-2 mm in thickness. The magnetic sus- 
ceptibility of a small piece of each of these specimens 
was measured using a standard Quantum Design SQUID 
magnetometer at Kyoto University. 

/iSR measurements were performed at TRIUMF, the 
Canadian National Accelerator Laboratory in Vancouver. 
Polarized positive muons were implanted one- by-one into 
the specimens mounted in a gas-flow He cryostat for mea- 
surements above T = 2.0 K at the M20 or Ml 5 channels, 
and in a dilution-refrgerator cryostat for those between 
15 mK and 10 K at the M15 channel. The time evo- 
lution A(t) of the muon spin polarization was obtained 
from the time histograms F(t) and B(t) of the forward 
and backward counters as 

A(t) = A G(t) = [F(t) - B(t)}/[F(t) + B(t)}, (1) 

where G(t) represents the relaxation function defined 
with G(0) = 1. Details of the /iSR methods can be found 
in refs4ii 



III. EXPERIMENTAL RESULTS: /iSR 
A. Zero-field /iSR spectra: ground state 

Due to its superb sensitivity to static magnetic order in 
spin systems with random / dilute / very small ordered 
moments, fiSR provides a very stringent test to verify the 
absence of active magnetism expected in spin-gap candi- 
date system s 26 ! 27 . Figure 3a shows the Zero-Field (ZF) 
/iSR time spectra obtained for (CuCl)LaNb2 07 which 
exhibit a very slow relaxation at T = 2 K and 15 mK. 
This depolarization can be decoupled by a small longitu- 
dinal field (LF) of 50 G. These features are expected for 
relaxation caused by static nuclear dipolar fields. To en- 
sure good heat conduction at 15 mK, we reproduced the 
results with another ceramic specimen containing 30% 
Au powder by weight. Thus, we confirmed the absence 
of static magnetic order in (Cu,Cl)LaNb2 07 down to 15 
mK. Statistical / systematic error of the measurement 
gives an upper limit of less than 2 % of the total volume 
Vm with static magnetism. 

In contrast, fast decay of the ZF-/iSR spectra was ob- 
served at T — 1.8 K in (Nbi-^Ta^) substitution systems 
with y > 0.3. With increasing the amplitude of the 
fast-relaxing component increases (Fig. 3a), indicating 
the existence of static magnetism with increasing volume 
fraction V M . In the CLAF system (CuBr)LaNb 2 7 , ZF- 
/iSR spectra A(t) exhibit long-lived precession below T/v 
(Fig. 3b) which indicates the existence of a well-defined 
local field at the muon site expected for homogenous 
long-range order. With decreasing Br composition x in 
the (Cli-^Br^) substitution, the internal field at T ~ 2 
K becomes increasingly inhomogeneous as shown by the 
damping of the oscillation in Fig. 3b. For x = 0.05, the 
inhomogeneous static local field results in a ZF-//SR line- 
shape often seen in dilute alloy spin-glass systems 2 ^. The 
change in the initial damping rate in Fig. 3b between x 
= 0.33 and 0.05 indicates that the spin structure / ori- 
entation and/or ordered moment size changes between 
these two concentrations. We confirmed a static origin 
of the observed fast relaxation in x = 0.05 via decoupling 
in LF at T = 15 mK. 



B. Zero-field /iSR spectra: temperature 
dependence 

Figure 4 shows the time spectra of ZF-/iSR observed in 
(CuBr)LaNb 2 7 and (CuCl)LaTa 2 7 at several different 
temperatures. Coherent oscillations are observed in the 
signal below the Neel temperature T^ = 32 K for the 
former and 7 K for the latter system. The spectra of 
(CuBr)LaNb 2 7 fit well to 

G z (t) = Ai[cos(cjt)] exp(-Ait)] 
+A 2 [exp(-A 2 t)] + A 3 [exp(-A 3 t] (2) 
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FIG. 3: (color-online) (a) Zero-field (ZF) /iSR time spec- 
tra in (CuCl)LaNb207 demonstrating absence of static mag- 
netism, and in (CuCl)La(Nbi- y Ta y )207 showing magnetic 
order in a partial volume fraction, (b) ZF-/iSR spectra in 
(CuCli- ;E Br ;E )LaNb2 07 at low temperatures. 
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FIG. 4: (color-online) Zero-field (ZF) /xSR time spectra in (a) 
(CuBr)LaNb207 and (b) (CuCl)LaTa2 07 showing the onset 
of long-range antiferromagnetic order. The solid lines repre- 
sent eq. 2 in (a) and eq. 3 in (b). 



where the first term represents the oscillating component, 
the second term represents a component showing a fast 
damping with A2 ~ 00 and the third term represents a 
persisting slowly relaxing signal with A3 << uj. The 
spectra of (CuCl)LaTa2 07 exhibit faster depolarization 
of the oscillation and fit well to a Bessel function term 
plus a nearly constant term, 

G z (t) = sm(cut) exp(-Ait)] 

+A 3 [exp(-A 3 £)]. (3) 

Previously, the Bessel function line shape was found in 
ZF-//SR spectra from systems having incommensurate 
spin-density- wave 29 or stripe spin modulation 4 . Neu- 
tron scattering measurements of (CuCl)LaTa207, how- 
ever, found a commensurate CLAF stated. Thus the 
present results may simply be due to highly-disordered 
short-range AF correlations. 

In Fig. 5, we show the observed frequency v = u/2tt 
for the (CI, Br) system with x = 1.0, 0.67, 0.33 and 



(CuCl)LaTa207. The spectra from the latter three sys- 
tems were fit to Bessel functions, in view of significantly 
better fits as compared to the cosine function, although 
none of these systems exhibit a clear indication of in- 
commensurate spin correlations in neutron scattering. 
With decreasing Br composition x, Tn decreases, but 
the frequency v(T — > 0) remains nearly unchanged at 
0.33 < x < 1. This indicates that the size of the or- 
dered Cu moment does not depend on x. The frequency 
v(T — > 0) for (CuCl)LaTa207 is significantly different 
from that of the other systems in Fig. 5, despite the fact 
that the same ordered Cu moment size ~ 0.6/i# was re- 
ported by neutron scattering studie s 19 ! 20 . The lower fre- 
quency / internal field was also found in all the (Nb,Ta) 
systems with 0.3 < y < 1 (as discussed later) as well 
as in the (CI, Br) system with x = 0.05 (see Fig. 3b). 
Interestingly, these systems with the lower internal field 
all have a Neel temperature of T/v ^ 7 K. These obser- 
vations suggest that compounds with sufficiently reduced 
T/v have a short-ranged spin structure and a possibly dif- 
ferent direction of the ordered Cu moments compared to 
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FIG. 5: (color-online) Muon spin precession frequency v ob- 
served in Cu(Cli-xBr x )LaNb 2 7 with x = 1.0, 0.67 and 0.33, 
and in (CuCl)LaTa 2 7 



systems with higher T/y. 



C. Results in Weak Transverse Field: volume 
fraction of the magnetically ordered region 

To determine the volume fractions of regions with and 
without static magnetic order, /iSR measurements in 



weak transverse field (WTF) are quite useful, as shown 
in reS. The persistent oscillation amplitude in WTF re- 
flects muons landing in a paramagnetic or non-magnetic 
environment. Figures 6a and b show the precessing am- 
plitudes in (Nbi-^Ta^) and (Cli-^Br^) systems mostly 
in WTF = 100 G. Some of the data were obtained with 
WTF ~ 50 and 30 G due to limitation of available spec- 
trometers. A sharp onset of the CLAF ordered state 
is seen for the Br substitutions with x = 0.2 - 1 (Fig. 
6b) at temperatures consistent with the magnetic sus- 
ceptibility results reported in paper III 19 . The (Nb,Ta) 
systems show static magnetism below a common onset 
temperature Tn ~ 6-7 K, with the paramagnetic volume 
fraction decreasing gradually with decreasing tempera- 
ture towards the partial fraction dependent on y. This 
confirms the phase separation observed in ZF-/iSR. In 
both Figs. 6a and b, about 25% of muons remain in a 
paramagnetic environment even when static magnetism 
is established in the full volume fraction as in the CLAF 
state of (CuBr)LaNb207^. These muons are likely oc- 
cupying sites where local fields from antiferromagnetic 
Cu spins nearly cancel for symmetry reasons. Judging 
from the amplitude from the remaining 75% of muons, 
we conclude that systems with x > 0.2 order in a full 
volume fraction Vm = 1.0, while those with x = 0.05 and 
0.3 < y < 0.9 undergo phase separation between ordered 
and spin-gapped volumes. The spin-gap volume prevails 
to nearly the full fraction for y < 0.3. 
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FIG. 6: (color) (a) and (b): Amplitude of persisting muon spin precession in a weak transverse field (WTF) of 100 G in (a) 
(CuCl)La(Nb,Ta)207 and (b) Cu(Cl,Br)LaNb207. This represents the fraction of muons landing in para- or non-magnetic 
environment. In each specimen, about 25% of the amplitude persists at any temperature and composition as denoted by the 
broken line, which presumably comes from muons at sites where the local fields from Cu moments cancel via symmetry reasons. 



D. ZF-/xSR spectra in phase-separated systems 



Static magnetism in the phase-separated region can in the (Nb,Ta) systems. The similarities of these two 
be elucidated by ZF-/iSR spectra in Figs. 7a and b for figures demonstrates that near the phase boundary be- 
quantum (Fig. 7a) and thermal (Fig. 7b) evolution tween CLAF and spin-gapped states, both quantum and 
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FIG. 7: (color) System (a) and temperature (b) dependence of the /iSR time spectra in zero field in (CuCl)La(Nb,Ta)207, 
which exhibit signals of fast- and slow-decay components, with composition- and T-dependent amplitudes, (c) The exponential 
muon spin relaxation rate A of the fast decay component in (a) and (b). 



thermal transitions involve phase separation with gradual 
change of Vm as a function of temperature T and com- 
position y. These spectra can be decomposed into two 
components with fast and slow relaxation, respectively, 
having T- and ^/-dependent amplitude ratios. The expo- 
nential relaxation rate A of the fast component is nearly 
independent of T and as shown in Fig. 7c, which im- 
plies that the ordered regions in different T and y share 
common microscopic spin arrangements of Cu moments. 
The common decay rate, indicative of the same spin con- 
figuration, is also found in the x = 0.05 (Br, CI) substitu- 
tion system (Fig. 2c). The present experiment, however, 
does not allow us to estimate the typical size of the or- 
dered regions. 



IV. EXPERIMENTAL RESULTS: LOW-FIELD 
MAGNETIC SUSCEPTIBILITY 

Motivated by the ZF-/iSR line shapes in Fig. 7 which 
resemble those expected in spin-glass systems 28 , we per- 
formed measurements of dc-magnetization M in a weak 
field of 100 G with both Field Cooling (FC) and Zero- 
Field Cooling (ZFC). As shown in Fig. 8a, marked depar- 
ture of Mpc from Mzfc sets in at T — 6 K, coinciding 
with the onset of static magnetism in the (Nbi-^Ta^) 
systems. History dependence of M was also found in 
(Cli-^Br^) systems with x < 0.66, with a common onset 
temperature of 6 K, well below XJv for x = 0.66, 0.33 and 
0.2 (Fig. 8b). The magnitudes of (M fc -M Z fc) at T = 
2 K for all of these systems correspond to a very small 
ferromagnetic polarization < 10 -3 /i# per Cu (insets of 
Fig. 8a and b), nearly independent of the field for the 
FC measurements in 100 - 800 G in the x = 0.05 system 
(inset of Fig. 8b). 

The irreversible magnetization M^ rr = (Mfc — Mzfc) 
in (CuCl)La(Nbi_ 2y Ta 2/ ) 2 7 at T = 2 K (Fig. 8a) roughly 
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FIG. 8: (color) (a) and (b): Magnetic susceptibil- 
ity for both Field Cooling (FC) and Zero-Field Cooling 
(ZFC) in 100 G obtained in (a) (CuCl)La(Nb,Ta) 2 7 and 
(b)Cu(Cl,Br)LaNb207- The inset figures show the irre- 
versible magnetization Mi rr = Mfc — Mzfc at T — 2 K. 



scales with the volume fraction Vm of the magnetically- 
ordered region (Fig. 6a). M irr and Vm share the same 
onset temperature. These features suggest that the ob- 
served signal likely comes from the bulk volume of the 
region with static magnetism, rather than from dilute 
ferromagnetic impurities. The very small net ferromag- 
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netic polarization (inset in Fig. 8a), together with the 
static Cu moment size of ~ 0.6 indicate mostly an- 
tiferromagnetic or random spin configurations, with a 
very small ferrimagnetic / canted component. Without 
distinguishing between these, we term the static mag- 
netism in the (Nb,Ta) system as a "glassy antiferro- 
magnetic" (GAF) state. The relatively large M^ rr ob- 
served in the y=l pure Ta material rules out the no- 
tion that the irreversibility requires randomness and/or 
a non-stoichiometic solid solution. Small, yet somewhat 
surprising non-zero values of M irr in the y = and 
0.2 systems may be due to increasing ferrimagnetic / 
canted contributions, possibly from remaining static re- 
gions with a volume fraction < 2%. 



V. PHASE DIAGRAM 

Summarizing these findings, we present a phase dia- 
gram in Fig. 9 as functions of the substitution concentra- 
tions x and y. The phase-separated region with a partial 
volume fraction is illustrated by a striped pattern, while 
regions with a full volume fraction are indicated by solie 
colors. The first-order thermal transition is shown by 
the broken line, while the solid line indicates the tran- 
sitions which are likely second-order. The internal field 
observed in the CLAF state (pink) is significantly larger 
than that in the glassy antiferromangetic (GAF) state 
(green), suggesting that these states are distinct. Further 
characterization is required to clarify the coexistence of 
these two states in the region of .05 < x < 1. The pur- 
ple and orange arrows indicate, respectively, the history- 
dependent (HD) region where Mpc 7^ Mzfc and the 
region with phase-separated static magnetism in a par- 
tial volume fraction (PV). 



VI. DISCUSSION 

The present /iSR results have demonstrated that 
(CuCl)LaNb2 07 indeed posseses a non-magnetic ground 
state, consistent with spin-gap formation, over the full 
volume fraction. /iSR has a superb sensitivity to static 
magnetism, as even small nuclear dipolar fields can easily 
be detected. Previous /iSR measurements on some sys- 
tems widely conceived to have spin-gap ground states, 
such as SrCu 2 (B0 3 )2 3 °, CaV 2 5 and CaV 4 9 31 , re- 
sulted in the detection of muon spin relaxation persist- 
ing to low temperatures. Contrary to these cases, the 
present study established the absence of any detectable 
static and dynamic magnetism from Cu moments, and 
demonstrated that the gapped state in (CuCl)LaNb2 07 
is really robust. 

Our results also revealed phase separation be- 
tween volumes with and without static magnetism in 
(CuCl)La(Nbi_ y Ta 2y )207 for the Ta concentration y > 
0.3. In the companion paper II, Kitada et al. 20 report 
neutron scattering results of the (Nb,Ta) systems, where 
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FIG. 9: (color) Phase diagram of Cu(Cl,Br)La(Nb,Ta) 2 7 ob- 
tained in the present study, with spin-gap (SG), glassy anti- 
ferromagnetic (GAF) and collinear antiferromagnetic (CLAF) 
states. The striped region indicates phase separation. The 
broken and solid lines denote, respectively, first- and second- 
order thermal transitions. The arrows attached to the hori- 
zontal axis show the region with history dependence (HD) of 
the magnetic susceptibility, and the region where static mag- 
netism exists in a partial volume fraction (PV). 



the volume fraction was estimated from the intensities 
of the magnetic Bragg peaks. Kitada et al. also decom- 
posed the response of high-field magnetization into the 
"gapped" and "ungapped" signals, using the results of 
the y = 0.4 compound as the reference for the former 
and y = 1 for the latter, and estimated the volume frac- 
tion of the ungapped signal. The volume fraction values 
derived from the neutron and magnetization results in 
these procedures agree well with Vm from /iSR. 

Both neutron and magnetization studies detect mag- 
netism as a volume-integrated quantity, and they can- 
not distinguish a small volume with a large individ- 
ual moment versus a large volume with a small mo- 
ment. In contrast, /iSR and NMR produce distinguish- 
able signals from magnetically-ordered and paramagnetic 
regions, with the amplitudes proportional to correspond- 
ing volumes: the volume information is decoupled from 
that of the moment size. Consequently, the real- space 
probes /iSR and NMR have genuine advantages over neu- 
tron and magnetization in the determination of ordered 
volume fractions. In the present (Nb,Ta) systems, the 
size of the ordered moment does not depend on the Ta 
concentration y, as demonstrated in Fig. 7. Due to this 
feature, the neutron and magnetization results for the 
volume fraction agreed well with the /iSR results. 

In paper III, Tsujimoto et al. 19 find signatures in 
the magnetic susceptibility for the (CI, Br) systems from 
which they derived the Neel temperature quantitatively 
consistent with the present /iSR results. Since the posi- 
tive muon is a charged probe, there remains some suspi- 
cion that the /iSR results may be different from those of 
the bulk system due to possible perturbation caused by 
the presence of the muon. The good agreements of the 
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FIG. 10: (color) A schematic view of free energy profile as 
a function of magnetic order parameter m. The curves (a) 
and (b) represent the case of standard second-order phase 
transitions, above and below T c , respectively, while (d) and 
(c) for the first-order transitions with (d) above T c and (c) at 



volume fraction and Tjy with the neutron and suscep- 
tibility / magnetization results confirm that the muons 
are indeed probing bulk properties in the present sys- 
tems. The /iSR results in Fig. 3(b) indicate increas- 
ingly short-ranged spin correlations with decreasing Br 
concentration x. This feature was not detected by neu- 
trons, due presumably to limited instrumental resolution 
and the weak signal expected from the magnetic Bragg 
reflections for powder specimens. The spin correlation 
length remains to be determined by neutron studies in 
the future using single crystal specimens. 

History dependence (HD) of the magnetization M has 
also been observed in the Kagome lattice antiferromag- 
nets SrCr ;E Gai2-xOi9 (SCGO) 32 , Cr-jarosite^ 3 - and many 
other GFSS, but without a jump in Mpc at the onset 
temperature. The Mpc jump observed in the present 
J1-J2 systems may be related to the involvement of the 
ferromagnetic exchange interaction Ji, which does not 
exist in other GFSS based exclusively on antiferromag- 
netic couplings. The partial volume fraction (PV) of the 
magnetically-ordered region adjacent to the spin-gap re- 
gion has also been observed recently by /iSR studies in 
the Kagome lattice Herbertsmithite system 26 . Figure 
10 shows a schematic view of free energy variation. A 
free-energy profile with multiple minima at the order- 
parameter (m) values of zero and a finite value in Fig. 
10 (lines c and d) can explain the origin of a first-order 
transition which often results in phase separation. His- 
tory dependence can be caused by multiple free energy 
minima, as generally seen in spin glasses with compli- 
cated free-energy landscapes with many minima. The 
present results show, however, that the regions for these 
two phenomena (purple and orange lines in Fig. 9) do 
not completely overlap. 



When we look into previous /iSR and susceptibility 
results in frustrated spin systems, we notice three pat- 
terns: (a) Spin freezing at T g , associated with crit- 
ical slowing down (maxima of 1/T\ of /iSR), history 
dependence of M below T p , and disappearance of dy- 
namics at T — » 0: canonical spin glasses AuFe and 
CuMn 28 . (b) Slowing down of spin fluctuations to- 
wards T g , history dependence of M below T g , fol- 
lowed by persistent dynamic effects at T — > 0, M 
remaining finite at T — > 0: Kagome lattice systems 
SCGO SrCr 8 Ga 4 0i 9 34 , Cr-jarosite KCr 3 (OH) 6 (S0 4 ) 2 33 , 
and volborthite (CuaZni-^V^CMOEb) 2H 2 35 , and a 
body-centered-tetragonal (BCT) system CePt2Sn2 36 . (c) 
Fully gapped state without static or dynamic magnetism 
at T — > with clear reduction of M at T — > suggest- 
ing spin-gap formation in the unperturbed system, and 
appearance of phase-separated static magnetism when 
composition is varied: CuClLaNb207 and the Kagome 
system Herbertsmithite Zn^Cu^^OH^Cb^. 

We also note that many of the "spin-gap candidate" 
systems, such as SrCu 2 (B03)2^, CaV^Os and CaV^Og^ 
and a doped Haldane gap system (Y,Ca)2BaNi05 3 - ex- 
hibit responses of the pattern (b), while the results of 
cuprate 4 ^ and FeAs superconductors 3 -^ 3 -^ correspond to 
the pattern (c). These observations indicate that near the 
fully spin-gapped (non-magnetic and/or superconduct- 
ing) state, systems have two choices: either to phase sep- 
arate and exhibit magnetic order in a partial volume frac- 
tion (pattern (c)) or to become a strange "spin-liquid" 
which shows persistent dynamic spin responses at T — > 
in the whole volume (pattern (b)). Many of the /iSR 
spectra observed in pattern (b) exhibit strange Gaussian 
line shapes which are very hard to decouple by longitu- 
dinal fields. Further studies are required to clarify which 
parameters are essential for a given system to follow pat- 
tern (b) or (c), and to characterize more details of spin 
dynamics and line shapes for pattern (b). 

Recent fiSR results in MnSi and (Sr,Ca)RuOari re- 
vealed phase separation at QPTs, very similar to the 
present results. Some history dependence was also no- 
ticed in MnSi 40 near the phase boundary where static 
magnetism disappears. Slow and presumably quantum 
spin fluctuations persisting at very low temperatures 
have been observed in MnSi 2,3 near the quantum phase 
boundary between the helically ordered and paramag- 
netic states. Phase separation in HTSC cuprates was 
observed in the QPT between spin-charge stripe and su- 
perconducting states 4 ^^. The magnetic resonance mode 
in HTSC is a quantum slow spin fluctuation and a soft- 
mode related to competing states across the quantum 
phase boundar y 41 ! 42 . Although a conclusive argument 
requires further accumulation of results on history depen- 
dence and inelastic soft modes in various systems, these 
observations hint that first order transitions, phase sepa- 
ration, history dependence, slow soft-mode spin fluctua- 
tions, and multiple free-energy minima may be common 
phenomena generic to QPTs in systems both with and 
without geometrical frustration. Comprehensive theoret- 
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ical studies on this aspect may reveal further fascinating 
features. 
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